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ABSTRACT: The mean area of two-dimensional unpressurised vesicles, or self-avoiding 
loops of fixed length A, behaves for large A as Aq A 3 / 2 , while their mean square radius of 
gyration behaves as Rq A 3 / 2 . The amplitude ratio Aq/Rq is computed exactly and found 
to equal 4tv/5. The physics of the pressurised case, both in the inflated and collapsed 
phases, may be usefully related to that of a complex 0(n) field theory coupled to a U(l) 
gauge field, in the limit n — > 0. 



Self-avoiding loops, or lattice polygons, have been studied extensively as models for 
planar vesicles. In a pioneering paper, Leibler, Singh and Fisher 1 (LSF) considered the 
statistics of the area and shape of such loops, subject to an osmotic pressure difference p, 
and controlled by a rigidity parameter k. While much of the interesting physics arises as a 
result of the variation of this latter quantity, these authors also observed interesting scaling 
behaviour as a function of p, when k = 0. Specifically, they found, on the basis of Monte 
Carlo studies and exact enumerations, that for an ensemble in which the total length, or 
mass, N of the loop is fixed, the mean area and squared radius of gyration behave as 

(A) N ~ A N 2 ^Y(x), (R 2 G ) N ~ R$N 2v X(x) (1) 

where x = pN^ v ', and v = | is the usual self-avoiding walk exponent in two dimensions. 2 
LSF argued that ip = 2, and conjectured vj± = is, a result which was derived indirectly 
by Duplantier 3 on the basis of Coulomb gas arguments. These results were confirmed 
and extended to measures of the shape dependence by Camacho and Fisher, 4 and the 
lattice enumeration studies were carried to higher orders in papers by Conway, Enting, 
Fisher, Guttmann and Whittington. 5 ' 6 ' 7,8 One result of these studies 4 ' 6 ' 9 was the apparent 
universality of the ratio Aq/Rq. 

This Letter describes an analytic approach to this problem. Generalising the well- 
known correspondence of de Gennes, 10 the problem of self-avoiding loops at fixed p and 
fixed monomer fugacity u (conjugate to N) is shown to be equivalent to a complex 0(n) 
spin model coupled to a U(l) gauge field. This field-theoretic formulation of the problem 
immediately establishes the scaling forms in (1). Moreover, using the methods of two- 
dimensional conformal field theory and the Coulomb gas mappings of Nienhuis, 11 the ratio 
Aq/Rq may be computed exactly to be 4iv/5. 

For the sake of definiteness, consider oriented self-avoiding polygons on a honeycomb 
lattice. The correspondence to a (complex) 0(n) spin model 10 ' 11 is as follows. Suppose 
that s a (r) label the components (a = 1, . . . , n) of a complex- valued spin at the site r. These 
spins are normalised so that Tr s*(r)sj,(r / ) = S a }jd rr i and Tr s a (r)si)(r') = Tr s^(r)s^(r') = 
0. Then the partition function Z = Tr Ylrr 1 0- + u E« s a( r ')' s tt( r ) +cc), gives, in the 
limit n — > 0, the generating function for the number of (unoriented) self-avoiding loops 
per site: Z = l + 2nftf^2 N pj^u N + (9(n 2 ), where Af is the total number of sites. Up to this 
point, the complex 0(n) model is completely equivalent to the usual real 0(2n) model. 
Consider now a unit current flowing along each link of a polygon in the direction of 
its orientation. 12 ' 9 ' 13 An explicit expression for this current, when inserted into correlation 
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functions of the lattice 0(n) model, is J^(r, r') = u{r' — r)^ J2 a i s a( r/ ) s a(r) — ex.), where 
the lattice spacing has been taken equal to unity. In the continuum spin version of the 
0(n) model, in which the spins are replaced by a field & a (r), Jfi is just the £7(1) current, 
(l/2z) J2 a (^ad^a — $ a ^^a)> whose space integral generates the global £7(1) symmetry 
$ a (r) — > e m $ a (r). The area of a given loop is now given by 

A = - l -J\ n - ri|*(r - rj) J (r) J (r')d 2 rd 2 r' (2) 

introducing Cartesian co-ordinates r = (ro,ri). This expression is valid for any non-self- 
intersecting loop (but not, in general, for loops which do self-intersect, since it weights 
different regions by the modulus of the winding number of the loop around them.) It 
differs from that for the signed area, which is proportional to e^ v J r jJj J l/ (r)d 2 r. This 
latter quantity was used in Ref. 14 to study ordinary intersecting loops as a model of 
vesicles. However, in the absence of any pressure difference its mean, after averaging over 
orientations of the loop, vanishes identically, and it is therefore not a suitable measure of 
the area. 

In fact, (2) is readily recognised as the expectation value of the polygon regarded as a 
Wilson loop in a £7(1) gauge theory: 

A = j (A^r)A v {r'))J^r)J v {r')d 2 rd 2 r' (3) 

in the gauge A\ = 0. It is straightforward to check that a similar result holds also in a 
covariant gauge, where (A^{r)A v {r')) = (l/2n){5 ^ v \iv\r — r'\ — r^r^/r 2 ^. This is impor- 
tant, since in this gauge the rotational invariance of the final result in the continuum limit 
is manifest. 

In the continuum limit near its critical point, the 0(n) lattice model corresponds to a 
field theory with 0(n) symmetry, which has been studied extensively in two dimensions. 15,16, 
This continuum theory will possess a conserved £7(1) current J^. It then follows from the 
above discussion that the continuum version of the generating function for self-avoiding 
loops weighted by is given by the n — > limit of an 0(n) theory, described by an 
action Sq, coupled to an Abelian gauge field: 

Z = Tr J VA^- So+ie I J » A ^ r -\ I F l^ r (4) 

Integrating out the gauge field yields the identification p = —e 2 /2. Several remarks may be 
made at this point. First, since A^ is dimensionless and J^, being a conserved Noether cur- 
rent, retains its canonical dimension of inverse length, it follows that p has renormalisation 
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group eigenvalue y p = 2. From this, and simple renormalisation group scaling arguments, 
follow the scaling laws in (1) with vj± = v and (p = 2. Second, the 'physical' region of 
the U(l) gauge theory, in which the gauge couping e is real and opposite charges attract, 
corresponds to a negative internal osmotic pressure difference. In that region, LSF find 
that, for large enough N, the loops collapse and behave like branched polymers. This is to 
be expected from the field theory, since in 1+1 dimensions a U(l) gauge field provides a 
confining potential so that the only asymptotic states are neutral. The world lines of these 
bound states correspond to the filaments of the branched polymer. The inflated phase, 
corresponding to p > 0, does not strictly make sense in the field theory, since the vacuum 
would become unstable to charge separation. As is well known, 17 the singularity in the 
free energy for e 2 < may be described in weak coupling by an instanton calculation. For 
p > it should be possible to neglect the self-avoiding constraint, so that the action Sq 
maybe replaced by that of n free complex scalar fields. In the first quantised picture, the 
instanton configuration corresponds to a particle-antiparticle pair being created at some 
imaginary time rg and annihilating at time 7q, their world lines describing a circle of radius 
R (fixed by extremising the total action) in euclidean space. This corresponds exactly to 
the physical picture of an inflated vesicle. 

Now return to the case p = 0, and the expression (2) for the mean area. It is easier 
to work in the axial gauge, although the same results are obtained in the covariant gauge. 
Averaging (2) over the ensemble of all self-avoiding loops, 

2nUY,PN{A) N u N = -\( |n - r[\5(r - r' )( J (r) J (r')>^V (5) 
N ZJ 

where J is now the U(l) current of the complex 0(n) theory, in the limit n — > 0. Thus 

i r°° 

nJ2PN(A) N u N = --a / ri(Jo(ri,0)Jo(0,0))dri (6) 

where ag is the area per site. In general, by current conservation and dimensional analysis, 
the correlation function ( J^(r) Jjy(0)), evaluated in the massive 0(n) field theory, has the 
form {d^d v — 5 ^d 2 ) f (m\r\) , where m is the mass and / is a dimensionless scaling function, 
whose large r asymptotic behaviour may be evaluated non-perturbatively using the form 
factor approach described in Ref. 13. However, at short distances, it becomes independent 
of m, and therefore has the form 12 k{n)(r^r v — \r^5^ v )/r^, determined, up to the con- 
stant k(n), by current conservation and rotational symmetry. Since the normalisation of 
this current is fixed by the requirement that its integral generate the U(l) symmetry, the 
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number k(n) is universal. Moreover, being a short-distance limit, it should be calculable 
within the conformal field theory corresponding to the massless complex 0(n) field theory. 
In field theory, k is called the chiral anomaly, since current conservation implies the exis- 
tence of a contact term, proportional to 5 2 (r), in the operator product Jz,(r) Jft(O) of the 
left- and and right-moving currents J^,_R = Jo^Jl- (For integral values of n > 1 the £7(1) 
symmetry would be embedded in a Kac- Moody algebra and k(n) would be proportional 
to the level number, but for n < 1 such a concept does not appear to make sense.) 

It then follows that the correlation function on the right hand side of (6) behaves like 
—k(n)/2(rf), so the integral appears to diverge logarithmically at short distances. In fact 
such a divergence must occur, since otherwise the integral would be dimensionless and 
therefore independent of the mass m, implying that the left hand side has no singularity 
as a function of u. In fact, if b is used as a short-distance cut-off, this divergence must 
be of the form ln(m&). Using the fact that the mass vanishes at the critical point u = u c 
according to m ~ (u c — u) u , it follows that the right hand side of (6) has the singular 
behaviour — ^aQk(n)u ha(u c — u) and hence that, as N — > oo, 

PN(A) N ~ \aa Q k'{V)vN- l u- N (7) 

where the lattice-dependent integer a appears 18 ' 9 because, on non-close packed lattices, 
the series is in fact in u a and therefore has a equivalent singularities on \u\ = u c - Defining 
the amplitude B by ~ BN^ v ^ l u^ N , it follows that BAq = ^aaQk f (0)u, which gives 
f fc'(O) for the square lattice (a = 2, uq = 1). 

The next step is to evaluate k(n) using Coulomb gas methods. 11 The mapping to the 
Coulomb gas proceeds in two stages. On the honeycomb lattice, the expansion of the 
partition function of the complex 0(n) model yields a sum over configurations of non- 
intersecting oriented loops, weighted by a factor of u for each link and n for each loop. 
This latter factor may be written in a local fashion by incorporating a factor e ±ix at each 
occupied vertex, depending on whether the oriented loop makes a turn through ±7r/3 at 
that point. Thus (anti-)clockwise loops accumulate factors of e~ 6 * x or e 6 ** respectively. 
After the summation over orientations, the appropriate factor of 2n per loop may be 
recovered by choosing n = cos 6%. This model may then be mapped onto a solid-on-solid 
(SOS) model by assigning heights <j)(r) (which are conventionally chosen to be integer 
multiples of n) to the sites of the dual lattice. Neighbouring heights on either side of an 
oriented bond differ by ±7r, otherwise they are equal. By convention, the higher side is on 
the right, looking along the oriented bond. This model is then supposed 11 to renormalise, 
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in the long wavelength limit, onto a Gaussian model with action Sq = (g/^ir) J"(<9</>) 2 <i 2 r, 
where g = 2 — 6x/ n - However, there is a caveat: in the SOS model the factors e ±%x lead to 
the result ( e - 12i ^MA) = 1; 

as may be seen by direct calculation in the fugacity expansion. 
In the Coulomb gas language, where 4> is interpreted as an electrostatic potential, this 
phenomenon corresponds to a total electric charge 12x/7r on the boundary, which preserves 
overall neutrality. Thus all non-zero correlation functions must correspond to a total charge 
— 12x/7r in the interior. 

How should the U(l) current of the complex O(n) model be represented in the SOS 
model? The naive candidate is simply = (1/ 7r)e^f/A^(^), where A;/ denotes a lattice 

difference between sites of the dual lattice. This current is conserved and has the property 
of taking the values ±1 as required. But it is nevertheless incorrect, since clockwise and 
anticlockwise loops are counted with different phase factors, resulting in a net clockwise 
current of e * — e~ * = 2i sin 6% around each loop. In addition, it may be seen that the 
correlation function (J^ (r)J„ (r')) receives contributions when the links r and r' are 
on different loops, a feature which is absent in the O(n) model. This correlation function 
also suffers from having net charge zero, so the charge on the boundary is not cancelled. 

However, there is another conserved current J^ 0S = Ae^Aj, ^e -12 *^/ 71 ^ , where the 
constant A is to be fixed. This has the property that its expectation value around a given 
loop vanishes, as required, on summing over both orientations, since this is proportional 
to (e~ 12 ** — l)e® lx + (e 12 ^ — l)e -6 ** = 0. Now consider the current-current correlation 
function, which is the expectation value of J^{r)J v {r'). For a given configuration in the 
loop gas, after summing over both orientations, this quantity takes the values 2nr^r' v if r 
and r' lie on the same loop (where f is a unit vector along the link r), and is zero if they 
lie on different loops. A suitable candidate for this in the SOS model is therefore 

(J% 0S {r)4° S {r')) = (A/7r)(e Ai7 A 7 e- 12 ^)/-e^A^(r / )) (8) 

since the quantity in (...) is zero when r and r are on different loops, for the same reason 
as above, and when they are on the same loop, it takes the value 7rf^f^(e _12 ^ — l)e^* — 
7r(e 12 *X — l)e~^* = —A-Kir^f' v sin 6%. Thus one should choose A = m/(2sin6x). It is 
somewhat curious that it is necessary to use different currents J and J in this expression, 
but such a result is enforced by the requirement that the total charge be — 12x/7r for the 
correlation function to be non-zero. 

It is now straightforward to evaluate (8) in the Gaussian model, replacing the lattice 
differences by derivatives, and using the fact that d^{r')e~ V2 ' ix ^l 1 ' ~ {-I2i x / gn){r - 
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- r / ) _2 e- 12i ^( r )/ 7r . The result is of the form expected, with 

. 12nx 2narccosn . . 

k(n) = — * = — (9) 

7r^sm6x 7rVl - n 2 (2n - arccosn) 

where < arccosn < it. For n = 1, the model describes a single species of charged boson 
with repulsive interactions, whose infrared behaviour is that of free fermions. In that case, 
one finds k = l/7r 2 , as may be checked independently. For n — > (9) gives k'(0) = 2/37T, 
so that BAq = aao/87r. 

In order to eliminate the lattice-dependent factors from this otherwise universal result, 
it may be combined with the relation BRq = 5cxao/327r 2 , which follows from a sum rule 1 ^ 
which is a consequence of Zamolodchikov's c-theorem 2 ^. This amplitude relation was first 
derived for the square lattice in Ref. 21, and generalised in Ref. 9. The main result given 
in the introduction then follows. The Table shows the comparison of the these predictions 
with results of lattice enumerations and a Monte Carlo simulation of a continuum model. 
The agreement is very satisfactory. 

It is interesting to note that random loops with self-intersection, for which (R 2 )]\r ~ 
N, correspond to a free 0(n) theory, for which the current- current correlation function 
behaves as r~ 2 lnr at short distances. 12 As a result, the mean area of such walks (weighted 
as discussed earlier) behaves as N\nN for large N. This is not surprising, since, as 
mentioned above, such loops are weighted by their winding number, whose average grows 
logarithmically. 22 In higher dimensions, (6) may be generalised to relate the generating 
function for the mean area of a loop projected onto a fixed plane to a similar integral over 
the current- current correlation function. Unfortunately, for d > 2 the singular behaviour 
of this integral does not come entirely from the short-distance behaviour, and therefore 
the whole scaling function f(mr) is required, rather than just the coefficient of its short- 
distance behaviour, which was calculated in an e-expansion by Miller. 12 However, this 
argument does imply that the mean projected area grows as N 2v ', in contradiction to the 
numerical findings (with rather short series) of Ref. 23. 

In this Letter I have shown how the theory of pressurised two-dimensional vesicles 
without rigidity may be given a field-theoretic basis. For zero pressure, this leads to an 
exact prediction for the universal amplitude ratio Aq/Rq. In the inflated phase p > 0, 
instanton techniques are applicable and should yield an asymptotic expansion for the 
crossover functions in (1) for large argument x. This work is currently in progress. The 
collapsed phase p < corresponds to confinement, and may provide an alternative field- 
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theoretic way of describing the so far only partially solved problem of two-dimensional 
branched polymers. 

I thank P. Fendley, M. E. Fisher, A. J. Guttmann and J. Miller for correspondence 
and discussions on this problem. 

Table. Most accurate existing estimates of the amplitudes B, Aq and Rq and their uni- 
versal ratios compared with the predictions of this work. 





Square 


Triangular 


Continuum 


Prediction 




2 


V3/2 






B 


0.5623 8 


0.2640 6 






A 


0.1416 5 


0.131 6 


0.314 ±3 4 




Rl 


0.05631 18 








BA /aa 


0.03981 


0.0399 




= 0.0397887 


Aq/RI 


2.515 




2.55 ±5 4 


3f = 2.51327 
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